arXiv: 1508.05150vl [gr-qc] 20 Aug 2015 


Saddle point inflation from f(R) theory 


Michal Artymowski 1 * * , Zygmunt Lalak 2 ! and Marek Lewicki 2 * 


1. Institute of Physics, Jagiellonian University 
Lojasiewicza 11, 30-348 Krakow, Poland 
2. Institute of Theoretical Physics, Faculty of Physics, University of Warsaw 
ul. Pasteura 5, 02-093 Warsaw, Poland 


Abstract 

We analyse several saddle point inflationary scenarios based on power-law f(R) models. We 
investigate inflation resulting from f{R) = R + a n M 2 0~ n ) R n -|_ a n +iM~ 2n R n+1 and f{R) = 
Y! n a n M 2 0~n) R n as well as l —> oo limit of the latter. In all cases we have found relation 
between a n coefficients and checked consistency with the PLANCK data as well as constraints 
coming from the stability of the models in question. Each of the models provides solutions which 
are both stable and consistent with PLANCK data, however only in parts of the parameter space 
where inflation starts on the plateau of the potential, some distance from the saddle. And thus 
all the correct solutions bear some resemblance to the Starobinsky model. 

1 Introduction 

Cosmic inflation HU2H3] is a theory of the early universe which predicts cosmic acceleration and 
generation of seeds of the large scale structure of the present universe. It solves problems of classi¬ 
cal cosmology and it is consistent with current experimental data 0] ■ The first theory of inflation 
was the Starobinsky model [5j, which is an f{R) theory [5] with R + R 2 /6M 2 Lagrangian density. 
In such a model the acceleration of space-time is generated by the gravitational interaction itself, 
without a need to introduce any new particles or fields. The embedding of Starobinsky inflation 
in no-scale SUGRA has been discussed in Ref. [7J. Recently the whole class of generalisations of 
the Starobinsky inflation have been discussed in the literature [mum iia Hug emi], also in 
the context of the higher order terms in Starobinsky Jordan frame potential HU EH ESI. 

The typical scale of inflation is set around the GUT scale, which is of the order of (10 16 GeU) 4 . 
Such a high scale of inflation seems to be a disadvantage of inflationary models. First of all 
inflationary physics is very far away from scales which can be measured in accelerators and other 
high-energy experiments. The other issue is, that high scale of inflation enables the production of 
super-heavy particles during the reheating [501 . Those particles could be in principle heavier than 
the inflaton itself, so particles like magnetic monopoles, which abundant existence is inconsistent 
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with observations, could be produced after inflation. Another argument, which supports low- 
scale inflation is the Lyth bound [21], which is the relation between variation of the inflaton 
during inflation in Planck units (denoted as A (f>) and tensor-to-scalar ratio r, namely 

< L1 > 

which for nearly scale-invariant power spectrum gives A <f> < M p for r < 0.002. Small Acf> seems 
to be preferable from the point of view of the naturalness principle, since M p is the cut-off scale of 
the theory. The value of r determines the scale of inflation, since V/r (where V is the potential of 
the inflaton) at the scale of inflation is set by the normalisation of CMB anisotropies. Therefore 
in order to obtain small r one needs a low-scale inflation, which may be provided by a potential 
with a saddle point. 

A separate issue related with f(R) inflation is related with loop corrections to the f(R) func¬ 
tion. In order to obtain quasi de Sitter evolution of space-time one needs a range of energies for 
which the R 2 M ~ 2 term dominates the Lagrangian density. This would require all higher order 
corrections (such as R 3 , R 4 etc.) to be suppressed by a mass scale much bigger than M. One 
naturally expects all higher order correction to GR to appear at the same energy scale if one 
wants to avoid the fine-tuning of coefficients of all higher order terms. From this perspective it 
would be better to generate inflation in f(R) theory without the Starobinsky plateau, which in 
principle could be obtained in the saddle point inflation. 

In what follows we use the convention 8nG = M p 2 = 1, where M p ~ 2 x 10 18 GeR is the 
reduced Planck mass. 

The outline of the paper is as follows. In Sec. [2] we give short introduction to f(R) and its 
description as a Brans-Dicke theory. In Sec. [3] we discuss three saddle point f(R) scenarios, 
namely: i) two higher order terms R n and R 11 ^ 1 , ii) at least 4 higher order terms with powers 
bigger than 2, iii) infinite number of higher order terms with finite sum at every energy scale. 
Finally we summarise in Sec. [4] 


2 Introduction to f(R) theory, inflation and primordial in¬ 
homogeneities 

The f(R ) theory is one of the simplest generalisations of general relativity (GR). It is based on 
Lagrangian density S = \ f d 4 ^/^gf(R) and it can be expressed using the so-called auxiliary 
field ip defined by ip = F(R) := In such a case the Jordan frame (JF) action is equal to 
S = f ^^g(ipR/2 — U(<p), where U = (RF — f)/ 2 is the JF potential. For F = 1 one recovers GR, 
so the GR vacuum of the JF potential is positioned at ip = 1. The same model can be expressed 
in the Einstein frame (EF), with the metric tensor defined by g pv = (pg pv . This is purely classical 
transformation of coordinates and results obtained in one frame are perfectly consistent with the 
ones from another frame^J The EF action is equal to S = f \J—g(R/2 + /2 — where 

R, (j> := ^8j2 logE and V := (RF-f)/(2F 2 ) are the EF Ricci scalar, field and potential respec¬ 
tively. The EF potential should have a minimum at the GR vacuum, which is positioned at (f> = 0. 


i 


Differences between Einstein and Jordan frame in loop quantum cosmology are described in Ref. 1221 . 
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In the EF the gravity obtains its canonical form and this is why the EF is usually used for the 
analysis of inflation and generation of primordial inhomogeneities. The cosmic inflation proceeds 
when both slow-roll parameters e and 77 are much smaller than unity. These parameters are, as 
usual given by 




( 2 . 1 ) 


where and are the first and the second derivative of the EF potential with respect to <p. 
During inflation e and 77 can be interpreted as deviation from the de Sitter solution for FRW 
universe. During each Hubble time the EF scalar field produces inhomogeneous modes with 
an amplitude of the order of the Hubble parameter. From them and from the scalar metric 
perturbations one constructs gauge invariant curvature perturbations, which are directly related 
to cosmic microwave background anisotropies. Their power spectrum V-r,, their spectral index 
n s and their tensor to scalar ratio are as follows 


'P'R. 


V 

247r 2 e ’ 


n s ~ 1 — 6 e + 277 , r ~ 16e. 


( 2 . 2 ) 


In the low scale inflation one obtains e <C |? 7 |, which for 77 < 0 gives 1 — n s ~ 2 |t 7 |. 


3 Saddle point inflation in power-law f(R) theory 

3.1 Saddle point with vanishing two derivatives 

As mentioned in the introduction, the loop corrections to the Starobinsky model are of the form 
L 2 ce n R n M 2 ( 1 ~ n ' ) , where M is a mass scale, which suppresses deviations from GR. Therefore, 
in order to obtain sufficiently long Starobinsky plateau one needs a broad range of energy scales 
on which R 2 dominates over all higher order corrections. This requires a fine-tuning of infinite 
number of a n coefficients. To avoid that we will consider an inflationary scenario in which 
different higher-order corrections can become relevant at the same energy scale, namely the 
saddle-point inflation from a power-law f(R ) theory. For general form of f(R) one obtains a 
saddle point of the Einstein frame potential for Va = Vaa = 0, which corresponds to 

RF =\f> RF’ = F , (3.1) 

where F = f and prime denotes the derivative with respect to the Ricci scalar. Let us assume 
the following form of f(R) 


f(R) - R + a + 


R n 




+ &n+l 


R 


n+1 


M 2n 


(3.2) 


where n > 2 is a given number. In such a case the saddle point appears for 

-1 

2 /'(n - 2)a n \ "- 1 f(n- 2)a 


R = R , = M 


n 


^n+l 


n 


(3.3) 


Equations above are a .2 independent, because any R 2 can satisfy Eq. (3.1). In order to keep 
R s and a n +i real we need to assume that a n > 0 and a n +i < 0. Then for sufficiently big 
R one finds F < 0 and the gravity becomes repulsive. This instability becomes an issue for 
R ~ M 2 ^ ( a n (n — 2)/n), which is typically of the same order of magnitude as R s . By 
redefining M we can always set one of a n to be any given constant. For negative a n one can 
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satisfy Eq. (3.1) for n < 2. Nevertheless the saddle point would he in the repulsive gravity 


regime, where F < 0. Thus in the following analysis n < 2 is excluded. Note that for non-zero 
value of «2 the value of M grows with < 22 . This comes from the fact that for 012 3> oi n one obtains 
inflationary plateau followed by the saddle point, due to growing value of R s with respect to < 22 - 
The «2 dependence of M is shown in Fig. |Tj Big <22 term means that the last 60 e-folds of 
inflation happen on the Starobinsky plateau, so one does not obtain significant deviations from 
the R 2 model. 


The Einstein frame potential around the saddle point (up to the maximal allowed value of 
<j>) for f(R) = R + «2 R 2 /M 2 + R 3 /M 4 + a^R A /M 6 has been shown in Fig [T] We have rescaled 
M to obtain a 3 = 1. The R 2 term in not necessary to obtain a saddle point, but we include it 
to combine the inflation on the Starobinsky plateau with the saddle point inflation. From Eq. 


(3.3) one finds the value of < 24 , normalisation of inhomogeneities gives M as a function of < 22 . 




Figure 1: Left panel: the Einstein frame potential for the (3.2) model around the saddle point for 
n = 3, «3 = 1 and different values of a 2 . The aq coefficient is set from Eq. (3.3) for n = 3. The 
maximal allowed value of cf is very close to the saddle point. Right panel: the scale of new physics 
M as a function of a. 2 - Dotted green line represent Starobinsky limit. 


3.2 Saddle point with vanishing k derivatives 


In general one can define the saddle point with first k derivatives vanishing, which was analysed 
in Ref. 233- In that case 1 — n s ~ N pf_i) when freeze-out of primordial inhomogeneities happens 
close to the saddle point. Thus, for sufficiently big k one can fit the Planck data. In our case 
all = 0 at the saddle point are equivalent to jpf = 0 for k > 2. The f(R) model from Eq. 


(3.2) cannot satisfy these equations, so in order to obtain a saddle point with vanishing higher 


order derivatives one needs to introduce more terms to f(R) function. Thus let us now consider 


f(R) -R + ar 

71=3 


R n 


M 2 (" _1 ) ’ 


(3.4) 


where l > 4 is an even natural number. Again, without any loss of generality one can choose <23 


to be any positive constant, so for simplicity we set <23 = 1. Then one can satisfy Eq. (3.1) and 
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Figure 2: Tensor to scalar ratio r and spectral index n s as a function of «2 for the (3.2) model with 
n = 3. One can fit the PLANCK data for «2 > 100, which means that the saddle point is preceded 
by the inflationary plateau. 


where 


/(") = 0 (for n = {3,4,..., l — 2} and any value of 02 ) and the saddle point appears at 
R = R S = y/pM 2 , 

The a n coefficients satisfy 

X„-1 2(1-3)! 


P=\l(l~ 1) ( 2 - 1 ) • 


- (-l) r 


(1 — n)\(n — 1)! 


for 


= {3 


(3.5) 


(3.6) 


Note that Eq. (3.5) and (3.6) are completely independent of 02 - Since on < 0 one obtains F < 0 


for sufficiently big R. Alike the model from Eq. (|3.2|) the biggest allowed value of R is slightly 

)'-l) 


bigger than R s . Using Eq. (3.4) and (3.6) one obtains 

f 


f(R) = R+^R 2 + R 


lM 2 flpR + 2 MV (( 




M 4 p — M 2 ^/pR 


(3.7) 


3.3 The l — > oo limit 

Numerical analysis shows that in order to obtain correct normalisation of primordial inhomo¬ 
geneities one needs M = Mil). Nevertheless for l —» oo one obtains M —>• M a (where M a ~ 10 -5 
for «2 =0), which implies R s —>■ oo for l — » oo. Hence for l 1 one cannot obtain inflation close 
to saddle point. For l —> oo one obtains 


f(R) = R 




(3.8) 


The «2 nray be again used to stabilise the GR vacuum at f = 0. The numerical results for 
N * = 60 are plotted in Fig. [6] and [tJ As expected, for a (» 1 values of M/y/a, r and n s obtain 
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Figure 3: Both panels present n s {ri) for the model from Eq. (3.2) for iV* = 50 and ct 2 = 0. If 
n is a natural number one cannot fit the Planck data due to too small n s . In all of those cases a 
significant contribution of the R 2 term is needed in order to obtain n s > 0.958. On the other hand 
for n — 2 < 10 -2 one obtain n s — 0.96. The case of n > 2 seems to be especially interesting since it 
allows to reconstruct Starobinsky results it the presence of higher order terms. 


the limit of the Starobinksy theory. As shown in Fig. [7] the potentials have two branches, which 
split at some <fi = </> m , where </> m is the minimal value of </>. The a 2 term in necessary in order 
to stabilise the GR vacuum. For a 2 = 0 one obtains two branches of potential which grow from 
(f> = 0. Both of them exist only for (f> > 0 with no minimum. While increasing the value of a 2 
the splitting of branches moves towards </> < 0 and the inflationary branch obtains minimum at 
(f> = 0. We investigated the stability of minimum from the perspective of classical evolution of 
the Einstein frame field. Namely, we considered the slow-roll initial conditions at f = </>* for 
different values of a 2 and checked whether the minimum is deep enough to stop the field before 
it would reach (f> m . We postpone the issue of quantum tunnelling to the anti - de Sitter vacuum 
for future work. 


4 Conclusions 

In this paper we considered several f(R) theories with saddle point in the Einstein frame poten¬ 
tial. All models consist of GR term R , Starobinsky term a 2 I? 2 and higher order terms which are 
the source of the saddle point. In subsection |3.1| we investigated two additional terms propor¬ 
tional to R n and R n+1 . We found analytical relation between their coefficients and R s , which is 
the value of the Ricci scalar at the saddle point. The potential becomes unstable for R slightly 
bigger than R s - the second branch of the auxiliary field equation = F(R) becomes physical, 
which leads to the second branch of potential and as a consequence to repulsive gravity. Signifi¬ 
cant contribution of the R 2 term extend the plateau before the saddle point and pushes away the 
instability from the inflationary region. For n > 3 it is impossible to obtain correct n s , however 
for n slightly bigger than 2 one can fit the PLANCK data. 

In subsection |3.2| we investigated The Einstein frame potential with zero value of the first 
l — 2 derivatives at the saddle point, where l > 6 is an even natural number. To obtain such 
a saddle point we considered f(R) = R + a 2 I ? 2 + Y(n=3 ot n M 2 P~ n lR n . We found analytical 
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Figure 4: Left panel: Numerical results for the model (3.1) for iV* = 50 and iV* = 60 (red and blue 
dots respectively). Right panel: Einstein frame potential for the model (3.1) for l = 6, l = 8, l = 10, 
l = 12 and l = 14 (orange, green, red, brown and blue lines respectively). The saddle point lies close 
to the right edge of the potential, beyond which one obtains a second branch of V, which leads to 
repulsive gravity. 


formulae for R s and for all a n coefficients, as well as the explicit value of f(R) after summation. 
Unfortunately the result is slightly disappointing, because the saddle point moves away from 
the scale of freeze-out of primordial inhomogeneities with growing l. Thus bringing us closer to 
the Starbinsky case as l gets bigger. We also obtained numerical results for n s , r and for the 
suppression scale M as a function of l. The final result strongly depends on N+, and therefore on 
the thermal history of the universe. One can fit the PLANCK data for l > 20 and IV* ~ 60 even 
for «2 = 0. Again, for R slightly bigger than R s one obtains an instability of potential, which 
for big l is orders of magnitude away from the freeze-out scale. 


In subsection 


3.3 


we considered the limit l oo, which resulted in f(R) = R(e _v/ ^ R / M ° + 
(y/2 + d‘ 2 )R/M )), which is basically Starobinsky model plus an exponentially suppressed cor¬ 
rection. In such a case the saddle point (and therefore the instability for R > R s ) moves to 
infinity and inflation happens far away from the saddle point. The oli term is necessary to create 
the meta-stable minimum of the Einstein frame potential. One can fit the PLANCK data for 
0.7 < «2 < 1.4 and a 2 ^ 34 
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Figure 5: Numerical results for the model (3.1) for N* = 50 and N* = 60 (red and blue dots 


respectively). All values of r obtained in this analysis are consistent with PLANCK, but n s fits the 
PLANCK data only for N+ ~ 60. 
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Figure 7: Left Panel: The Einstein frame potential as a function of the Ricci scalar. The GR 
minimum at R = 0 appears to he meta-stable, with a possibility of tunnelling to anti de Sitter 
vacuum. Right Panel: The Einstein frame potential V as a function of the Einstein frame field <ft 
for the model (3.8). Two branches of potential correspond to two solutions of ip = F(R). In order 
to avoid overshooting the minimum at R = 0 one requires a 2 > 0.7. 



Figure 8: The minimal value of 012 , which allows to avoid overshooting the meta-stable GR minimum. 
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